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自 适应 二 进 小 波 去 噪 法 
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摘 要 : 本 文通 过 引入 一 类 新 的 可 微 阔 值 函 数 ， 建 立 了 信和 号 的 能 量 同 它 的 二 进 小 波 变换 系数 之 间 的 关系 。 
利用 Stein 无 偏 风 险 估计 和 梯度 下 降 法 ， 给 出 了 迭代 求解 各 尺度 二 进 小 波 系数 最 优 阔 值 的 自 适 应 
去 噪 算法 。 数 值 实验 结 果 表 明 ， 同 现存 的 小 波 去 噪 技术 相 比 ， 该 算法 的 去 噪 结果 具有 更 高 的 重 构 


信 品 比 和 更 光滑 的 外 观 。 
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1 引言 


从 噪声 中 恢复 未 知 信号 一 直 是 信号 处 理 领 域 广泛 研究 的 重要 内 容 。Donoho 和 Jonestone 在 
小 波 变换 的 理论 基础 上 上， 创立 了 小 波 去 噪 法， 该 方法 能 有 效 地 剔除 信号 中 的 加 性 随机 噪声 ， 是 
信和 号 估计 的 一 个 简单 的 且 强 有 力 的 工具 ， 并 已 证 明 在 广泛 的 函数 空间 内 (如 Besov 和 Triebel 2 
间 等 )， 小 波 去 噪 的 效果 明显 优 于 其 它 经 典 的 去 品 方 法 (如 核 估计 法 、 样 条 平滑 法 和 Fourier 级 
数 展开 法 等 )， 并 且 在 高 斯 白 噪声 背景 下 ， 该 方法 在 均 方 误差 极 小 一 极 大 意义 下 是 几乎 最 优 的 
估计 器 H 3 。 

从 Donoho 和 Jonestone 的 开拓 性 论文 2 发 表 以 来 ， 人 们 在 小 波 去 噪 领域 已 做 了 大 量 的 研 
究 工 作 ， 其 中 大 多 数 工作 主要 集中 于 如 下 两 个 方面 : 

1.1 数据 自 适应 的 小 波 去 噪 法 

对 于 给 定 的 含 噪 信号， 小波 去 噪 的 效果 依赖 于 阔 值 函数 和 各 尺度 上 闪 值 的 选取 以 及 
噪声 的 实现 特别 ， 阐 值 的 选择 是 至 关 重 要 的 。 如 果 阔 值 取得 太 大 ， 如 象 通常 情况 下 那样 
取 上 = oV2log。 入 有 可 能 会 因 过 多 地 除去 信号 的 成 份 而 使 恢复 的 信号 出 现 严 重 的 失真 外， 而 阐 
值 太 小 又 难以 收 到 理想 的 去 噪 效果 。 因 此 ， 根 据 含 蝶 信号 样本 ， 分 尺 地 自 适 应 优化 阔 值 选择 是 
十 分 必要 的 。 在 这 类 去 噪 方法 中 ， 主 要 包括 基于 Stein 无 偏 风 险 估计 的 小 波 去 噪 算法 [8&6]、 自 适 
应 Bayes 小 波 去 噪 算法 [7,8] 和 块 阔 值 小 波 去 噪 算法 [19 等 。 

1.2 ”基于 宛 余 基 集 的 小 波 去 噪 法 

小 波 去 噪 的 性 能 主要 依靠 于 信和 号 在 给 定 小 波 基 下 稀疏 表示 的 程度 。 当 信号 包含 更 加 复杂 的 
结构 时 ， 经 过 二 进 制 伸缩 和 平移 生成 的 正 交 小 波 基 ， 通 常 不 能 较 好 地 匹配 信号 的 局 部 时 频 结 
构 ， 从 而 导致 信号 的 能 量 在 小 波 系数 域 中 扩散 。 在 这 种 情况 下 ， 常 常 导致 原来 低 于 阔 值 的 噪声 
小 波 系数 因 信和 号 成 份 的 加 入 变 得 高 于 阔 值 而 保留 下 来 ， 而 含有 重要 信息 的 信号 小 波 系 数 可 能 
低 于 阔 值 被 置 为 零 。 因 此 在 信号 重 构 时 ， 一 方面 会 使 噪声 过 多 地 留 在 重 构 信号 中 ， 另 一 方面 会 
因 信 和 号 能 量 的 损失 在 急剧 变化 部 分 产生 震荡 效应 ， 甚 至 会 使 重 构 信 号 发 生 严 重 变 形 。 因 此 ， 
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人 们 提出 了 小 波 包 去 噪 算法 、 平 移 不 变 的 小 波 去 噪 算法 或 非 抽取 小 波 去 品 算 法 上 3,1I44， 以 及 基 
于 Curvelet、Ridgelet 去 噪 算法 [15,16] 等 基于 宛 余 基 集 的 小 波 去 噪 法 。 

本 文 考虑 了 基于 二 进 小 波 变 换 的 阔 值 自 适应 选择 的 去 噪 问 题 ， 建 立 了 信号 的 能 量 同 它 的 二 
进 小 波 系数 能 量 的 关系 ， 给 出 了 一 类 新 的 可 微 的 阔 值 函数 ， 利 用 Stein 无 偏 风 险 估计 ， 采 用 梯 
度 法 对 各 个 尺度 上 的 阔 值 进行 自 适 应 地 优化 选择 ， 从 而 给 出 了 一 种 基于 二 进 小 波 变 换 的 自 适应 
去 噪 算法 。 数 值 实验 结果 证 实 了 该 算法 的 有 效 性 。 


2 二 进 小 波 变换 


本 节 首 先 介绍 二 进 小 波 变 换 的 定义 及 其 快速 分 解 与 重 构 算法 ， 然 后 给 出 信号 能 量 同 它 的 二 
进 小 波 变换 能 量 的 关系 ， 该 关系 对 下 一 节 自 适应 去 噪 算法 的 建立 是 至 关 重 要 的 。 
设 y(t) 为 小 波 函 数 ， 对 任意 信号 f(t)e LR), WES fo) 的 二 进 小 波 变换 定义 为 


Ws f(t) = F * Va (È), 
因此 ， 二 进 小 该 变换 是 对 连续 小 波 变换 的 尺度 参数 进行 二 进 制 离散 ， 是 一 种 平移 不 变 的 小 波 变 
换 。 为 了 确保 二 进 小 波 变换 的 重 构 公 式 成 立 ， 小 波 函 数 必须 满足 稳定 条 件 


4<》 Wio)| <B, 
. JEA 
其 中 常数 0 < 4 < B < oo。 这 时 ， 若 存在 函数 x(t) < L7(R), ACH Fourier 变换 满足 
》 b(2iw)x(2'w) = 1, 
jEZ 
则 称 x(t) A b(t) 的 一 个 重 构 小 波 。 这 时 信号 f(t) 可 从 它 的 二 进 小 波 变 换 获得 完备 和 稳健 的 重 构 
f(t) = > Wy f * Xz (0). 
jez 


一 个 离散 信号 二 进 小 波 变换 具有 快速 的 分 解 与 重 构 滤波 器 组 算法 ， 即 所 谓 的 多 和 孔 
算法 B719。 设 {hn]} 和 {gln]} 为 尺度 函数 p 和 小 波 函 数 力 对 应 的 共 轿 滤波 器 ， 令 aj[ln] = 
(F(t), 2-7/7 G9(2-F(t — n))), djin] = (FQ, 2-47 W(24(E 一 mn)))，fIm] 为 f(t) 的 离散 采样 ， 则 
对 任意 7 > 0， 离 散 信号 fm] 的 二 进 小 波 变换 快速 分 解 和 重 构 算法 如 下 。 


分 解 算法 : aj41[n] =a; *hj[n], djtiln] = a; *9,[n], (1) 
EHSA: ajin] = 5(az+1 * hyla] + diz * g;l]), (2) 


其 中 
ao 四 = fin), hjln] =h;[-n], 9,[n] = 9,[-nl, 


hi M gi 分 别 由 疡 和 9 在 每 对 相 邻 样本 间 插 2? — 1 个 零 得 到 。 
引 理 1 对 于 任意 尺度 sj > 0, Aljaz? + lld? = 2lle 1。 


第 6 期 KAE, fea: 自 适应 二 进 小 波 去 噪 法 971 


证 明 对 任意 尺度 j > 0, Wâl), dw), H) MG (w) SHH {asin}, {dsla}, (hinl) 和 
{g|n]} 的 离散 Fourier 变换 ， 则 有 


llagll? + [lag]? = Wall? + Udy |]? 


1 ft T 1 fr a 
af Bo Pao Pwt E f CaP 


1 us 
2 f las-s(w)Paw = 2 aj- IP. 


从 而 ， 引 理 得 证 。 
命题 1 对 于 任意 尺度 J > 0， 有 


J 
27 llaol? = llaz 要 十 > 27-7 45. 
j=1 
证 明 用 数学 归纳 法 证 明 。 
当 J 了 = 1 时 ， 由 引 理 知 ， 故 命题 成 立 。 假 设 当 J 二 上 时 命题 成 立 ， 即 
k 
2 laol? = llall? + > 2 lazh’, 
j=l 


则 当 J =k 十 1 时 ， 由 上 式 和 引 理 ， 有 


k k 
2 Jao? = 2(ljaxl? + 52-7 ldjll) = 2loxl? + D2 ay 
j=l j=l 


k k+1 
= llaky? + lidel? 十》 2**2-9 [dj lI? = larl? + 5° 2k? aN’. 


j=l j=l 


因此 ， 命 题 得 证 。 


3 自 适应 去 噪 算法 
3.1 一 类 新 的 可 微 阅 值 函数 
广泛 使 用 的 硬 立 值 函 数 是 不 连续 的 ， 从 而 阻碍 了 使 用 梯度 法 进行 自 适应 阐 值 选择 。 下 面 ， 
给 出 一 类 新 的 具有 连续 高 阶 导数 的 阔 值 函数 ， 它 们 可 以 任意 地 双 近 硬 阐 值 函数 。 令 
ge” l/r)? xz>t, 
6" (x,t) = 0, lz| <t, (3) 


ze—i/r(att)? r< 一 上 
其 中 ~ > 0。 从 上 述 8"(z,t) 的 定义 ， 立 即 有 lim 6"(a,t) = 6*(z,t)， 并 且 可 以 容易 地 求 得 


BO0"(z,t) 80"(z,t). O20"(z,t) _ 076"(z,t) 


Ot? Ox Orot OtOx 
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的 解析 表达 式 。 

E12 h T Br BY PB A eX o (a, t) (ME BR) Fl 4 (x, t) (RAR), DA De HB A R BLO" (a, t) ( 实 
线 )。 从 图 可 以 看 出 ， 随 着 7 的 增加 ，0"(zx,t) BOR BREATHE O” (x,t) Auk, ENARA 
相似 于 硬 阔 值 函数 的 作用 ， 同 时 ， 在 阔 值 的 附近 ， 它 亦 具 有 软 阔 值 函 数 的 作用 ， 它 可 使 得 小 波 
域 滤波 以 较为 平滑 的 方式 实现 ， 使 在 阔 值 的 附近 的 小 波 系数 缩减 后 更 加 接近 原 信号 的 小 波 系 
数 ， 从 而 可 抑制 由 直接 截断 可 能 引起 的 虚假 震荡 。 


图 1: BVA BS Ol (x,t). O4(x,t) FO" (a, t) HA 


3.2 ”去 噪 算法 
设 { 人 1 为 含 噪 的 观测 序列 ， 它 可 以 描述 为 如 下 模型 


yli) = fli] + nfi, O<i<N-1, (4) 


其 中 /加 为 未 知 确定 信号 的 采样 值 z 轩 是 独立 同 分 布 的 高 斯 白 噪声 ， 且 mlil ~ N(0,07). Æ 
记 


y = (yio, yi] = ,yIN -= 1])”, 


f= (flO), fl]... fIN-1])’, 
n= taia aN —1])7, 


则 式 (4) 可 重新 表示 为 
y =f+n, (5) 

于 是 ， 基 于 二 进 小 波 变 换 的 自 适应 去 噪 算法 可 描述 为 : 

(a) 二 进 小 波 分 解 

选择 适当 的 小 波 函 数 销 ， 并 且 确 定 分 解 层 数 .J(J < log N) {h[n]} Hl {g[n]} 分别 为 尺度 
函数 pg 和 小 波 函 数 几 对 应 的 共 轿 滤波 器 ， 对 观测 信号 y 使 用 式 (1) 作 离 散 二 进 小 波 分 解 (a。 = 
y)， 得 到 它 的 各 层 小 波 系数 dj = {d,,|0Sk<N}, 1<7J<J， 和 第 J 层 的 逼近 系数 av = 
{a,, |O<k< N}. 

(b) BEM BE A) R RURE 

使 用 (3) 式 定义 的 阔 值 函数 0"(z,t)， 分 别 对 各 层 小 波 系 数 进行 处 理 。 令 


dj = {djk =O" (dikti) |0<k<N}, 1<j<J, 
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HP t > 0 为 第 7 层 小 波 系数 的 阔 值 。 
下 面 ， 使 用 Stein 无 偏 估计 器 Q9 和 梯度 法 ， 给 出 各 层 小 波 系数 阔 值 的 自 适 应 选择 。 我 们 的 
HRE RERE = (fiz ,二 )T"， 使 它 极 小 化 重 构 均 方 误差 


R(t) = E( |f — f°), (6) 
其 中 dj = {df 10 < k < N] 为 信号 f 的 第 5 层 小 波 系数 ，1 < j <J, af = {af ||0<k<N} 


为 信号 f 的 第 了 层 的 逼近 系数 ,f= (FO), A0- , FIN -1)7, fin), 0<n< NIH, 1< 
JS J 和 a 使 用 重 构 公式 (2) 得 到 的 重 构 信 号 ，t = (tta ,ty)7。 根 据 第 2 节 中 的 命题 


式 (6) 可 重 写 为 


J 
R(t) = E (IE-t?) = B (D2 Af + 2 Jay — af?) 
2 ; ~ 2 2 
= 2029 E (lä; - afl?) + 2-78 (llas — af 7), 
j=1 


因此 ， 求 解 极 小 化 R(t) 的 问题 就 等 价 于 求解 
mE -afl), 1<jg7, 
F He, (djk) = 8 (dirt) — die = Gk -djr WU He, (-) 是 可 微 的 ， E EE E 
从 而 ， 据 Stein 的 原理 知 [19 
Rj(t;) = No? + 202VH;, (d;) + |H; (45) |[?, 
为 Elid; — df |)? 的 一 个 无 偏 估计 器 ， 其 中 


T N ƏH, (2) 
Hi, (dj) = (Hi, (d;o), Ht, (dja), » Ha; (dj,n-1)) V - H; (d;) = or 
k=1 


z=djk 
因此 ， 求 解 t = (6, i iT, AE 
t = arg min R(t), 
的 问题 ， 就 转化 为 求解 t* = (tt, th, --- ary dE 
= arg min Rott) j=1,2, AS. 


THe, MR (ts), 1 < j < J， 分 别 使 用 梯度 法 ， TS CSR AF tf, t3,--- ,她 ， 具 体 算法 描述 如 下 


O(N-1R,(t; 
H(t 1) = (8) | k 
其 中 oj (k) JERE K 
ƏR; (ts) N-1 8H, (2) ae 0? Hı (x) 
=2 Hi, (d; k) — +2 za 
Dt; a t; ( jk) at; T=d; k a 2 Orot; T=dj,k 


(c) 信和 号 的 恢复 
对 滤波 后 的 小 波 系数 du, 42,… ,Gy 和 逼近 系数 ay， 按 二 进 小 波 变换 的 重 构 公式 (2) 进行 重 
构 ， 从 而 得 到 去 噪 后 的 信号 f= (flo], A, ,FLV — 1p. 
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4 ”数值 实验 


使 用 本 文 提出 的 去 噪 算法 ， 我 们 对 Donoho 给 出 的 Blocks 信号 、Bumps 信号 、Heavesine 信 
号 和 Doppler 信 和 号 等 各 种 不 同时 频 结 构 的 信号 进行 了 数值 试验 ， 获 得 了 和 良好 的 去 噪 效果 。 
由 于 篇 幅 限 制 ， 在 下 面 的 实验 中 ， 仅 给 出 Bumps 信 和 号 的 去 噪 结果 。 该 信号 在 [01] 区 间 中 
采样 ， 采 样 点 数 为 N = 1024， 采 用 Daubechies 小 波 db8 作为 基本 小 波 ， 小 波 分 解 层 数 .J = 
6 用 来 比较 的 去 噪 算法 包括 Matalab 小 波 工具 箱 中 Sqtwolog 去 噪 算法 、Heursure 去 噪 算法 以 
及 Coifman 和 Donoho 提 出 的 平移 不 变 的 小 波 去 噪 算法 ， 它 们 均 使 用 硬 阐 值 函数 。 在 本 文 
提出 的 去 噪 算法 中 ， 使 用 由 式 (3) 定 义 的 阐 值 函数 94(z,t)， 在 阔 值 迭代 公式 (7) 中 ， 初 值 取 
Ai t;(0) = oV2logeN, 1 < j < J, BREKMA aj;(k) = 0.05, 1 < 7 < J k>0. EXMAM 
Bi i 


std({ fl] }) std({f[é]}) 
SN Rinput = 20logio ( 二 PE a k SN Routput = 20log10( aan J 
则 改善 信 品 比 为 SN Rimprove = SN Routput 一 SN Rinput， 其 中 std(.) 表示 序列 的 标准 方差 ，f 为 
了 的 去 噪 信号 。 

图 2 给 出 了 本 文 提出 算法 ( 实 线 )、 平 移 不 变 去 噪 算法 (点 线 )、Heursure 去 噪 算法 (虚线 ) 和 
Sqtwolog 去 噪 算法 (点 虚线 ) 的 信 噪 比 改善 曲线 。 考 虑 到 噪声 样本 的 随机 性 ， 对 于 每 个 给 定 的 输 
入 信 噪 比 ， 改 善信 噪 比 为 100 次 不 同 噪声 样本 下 的 改善 信 噪 比 的 均值 ， 从 图 2 可 以 看 出 ， 对 所 
有 输入 信 品 比 ， 本 文 提 出 算法 的 改善 信 噪 比 均 是 最 高 。 图 3 给 出 了 本 文 提 出 的 算法 在 输入 信 品 
比 SN Rinput = 20log106dB (15.56 dB) 时 的 迭代 收敛 曲线 ， 同 图 2 相同 ， 它 是 100 次 不 同 噪声 样 
本 下 的 平均 。 图 4 给 出 了 输入 信 噪 比 SNRinput = 20logjo6dB 时 不 同 算法 恢复 信号 的 效果 ， 从 
该 图 可 以 看 出 ，ADAPWDEN 算法 可 以 更 好 地 保存 原 信号 中 的 尖峰 和 突变 部 分 ， 同 时 也 更 好 
地 抑制 了 噪声 。 这 是 因为 在 ADAPWDEN 算法 中 ， 各 层 小 波 系数 的 阐 值 是 根据 极 小 化 重 构 均 
方 误差 自 适应 地 选取 的 ， 以 及 新 的 阔 值 函数 使 小 波 系 数 域 滤波 以 较为 平滑 的 方式 实现 ， 并 且 在 
几 值 附近 具有 相似 于 软 阔 值 函数 的 作用 ， 使 得 滤波 后 小 波 系数 的 取 值 更 加 接近 原 信 号 小 波 系数 
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的 取 值 。 
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小 波 去 噪 的 效果 主要 依靠 于 小 波 变换 、 阔 值 以 及 阔 值 函数 的 选择 。 在 这 篇 文章 中 ， 考 虑 了 
基于 二 进 小 波 变换 的 自 适应 阔 值 选择 的 去 噪 算法 。 首 先 ， 给 出 了 一 类 新 的 可 微 阐 值 函数 ， 这 类 
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阔 值 函数 可 以 任意 地 逼近 硬 阔 值 函 数 。 其 次 ， 建 立 了 信和 号 的 能 量 与 其 二 进 小 波 变换 系数 能 量 的 
之 间 关 系 ， 进 而 把 极 小 化 重 构 信 和 号 的 均 方 误差 问题 转化 为 分 别 极 小 化 各 层 小 波 系数 的 重 构 均 方 
误差 问题 。 最 后 ， 利 用 Stein 无 偏 估计 原理 ， 给 出 了 自 适应 求解 各 层 小 波 系数 最 佳 阔 值 的 梯度 
下 降 算 法 。 由 于 二 进 小 波 变换 的 平移 不 变性 和 宛 余 性 ， 以 及 各 尺度 小 波 系数 阔 值 选择 的 最 优 
性 ， 因 而 同 现 有 的 小 波 去 品 算 法相 比 ， 本 文 提出 的 去 噪 算法 不 仅 可 以 更 好 地 抑制 重 构 信号 中 的 
人 工 噪 声 ， 具 有 更 高 的 重 构 信 噪 比 和 更 光滑 的 外 观 ， 并 且 能 更 好 地 保持 原 信号 中 的 尖峰 和 突变 
成 分 ， 这 些 从 数值 实验 结果 中 得 到 了 验证 。 
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An Adaptive Denoising Algorithm Based on Dyadic 
Wavelet Transforms 


ZHANG Zhuo-sheng', REN Pin-yi? 


(1- School of Science, Xi’an Jiaotong University, Xi’an 710049; 
2- School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049) 


Abstract: A new class of thresholding functions is proposed. The energy conservation of the dyadic 
wavelet transform is proved. An adaptive denoising scheme based on the dyadic wavelet transform is 
introduced, which searches for the optimal threshold of each wavelet subband by minimizing a Stein 
unbiased risk estimation. Numerical experiments demonstrate that the proposed approach outperforms 
the classical wavelet-based denoising techniques. 
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